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The distribution of the capture zones of point islands nucleated in one dimension are considered using a 
retrospective view, providing an alternative perspective on why scaling occurs in this continuously evolving 
system. Distributional fixed point equations for the capture zones are derived through a mean field approx- 
imation for nearest neighbour size correlation. Solutions of these equations compare favourably to Monte 
Carlo data for a range of critical island size i = 0,1, 2, 3, and to theoretical models based on more traditional 
fragmentation theory approaches. 
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Scale invariance during the nucleation and growth of 
islands driven by monomer deposition is an intriguing 
phenomenon™. Island size distributions, and the distri- 
bution of capture zones which underlie the island growth 
rates, evolve towards scaling forms despite on-going nu- 
cleation of new islands with the concomitant disruption 
to the existing capture zones 2 -. The form of the scaling 
functions depends on the critical island size i, where i + 1 
is the smallest stable island size. A number of theoreti- 
cal approaches have been used to model this behaviour, 
ranging from mean field models which neglect the varia- 
tion in capture zone sizes^r— due to spatial arrangements 
of the islands, to those which attempt to include this 
information explicitly^—. All these approaches can be 
characterised as forward-looking in the sense that they 
are based on predicting how size distributions evolve as 
new islands nucleate. In this paper we shall present 
an alternative, retrospective, perspective where we ask 
how the capture zones present in the system came to 
be created. We focus on the case of point island nu- 
cleation in a one dimensional (1-D) system since this 
allows for a more complete analysis and comparisons 
with results from a more traditional fragmentation theory 
approac h 14 ! 15 . This new perspective provides interesting 
insight into why scaling occurs and compares well with 
simulation data. 

Island nucleation and growth is widely studied using 
Monte Carlo (MC) simulation. A point island approxi- 
mation is often used both for clarity and because it ap- 
proximates the growth of small, well-separated islands 4 . 
Here we employ a 1-D model 14 where monomers are de- 
posited at random onto an initially empty lattice at a 
deposition rate of F monolayers per unit time. The 
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monomers diffuse at rate D on the lattice, nucleating 
immobile point islands when i + 1 monomers coincide at 
a lattice site. Once nucleated, the islands grow by ab- 
sorbing any monomers that hit them, although they only 
occupy one lattice site. When sufficient islands have been 
nucleated, the most likely fate of a deposited monomer 
is to become absorbed by an existing island rather than 
being incorporated into a new island. It is in this aggre- 
gation regime of growth where scale invariance is found; 
note however that island nucleation continues still, albeit 
at a slow rate compared to monomer adsorption. 
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FIG. 1. The islands numbered Ji — J5 on the 
one dimensional substrate. The gaps between 
the islands are labelled gi, gs, g$ and 34, and 
the capture zones of islands I3, I5, I4 are la- 
belled C3, C5 and C4 respectively. 



Figure Q] shows some islands on the lattice, numbered 
according to their chronological age, along with their cap- 
ture zones C3, C4 and C5. Island I3 has the capture 
zone of size C3 = (g% + gs)/2, where g\ and (73 are the 
inter-island gaps to the left and right of ^3 respectively. 
C3 represents the average growth rate of I 3 , since any 
monomers deposited into C3 are more likely to diffuse to 
I 3 than its neighbours I± and I 5 . 

Referring to Figure [TJ let us ask how the inter-island 
gap 173 was created. It was formed by the nucleation of 
the youngest island in the picture, 1$, which occurred 
in the gap of size (33 + #5) between islands I3 and I a- 
Generalising, we will suppose that any randomly chosen 
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gap with size x (scaled to the average) in the system will 
have arisen by the fragmentation of a larger gap formed 
by combining it with a neighbouring gap of size y. In gen- 
eral we do not have the benefit of the chronological ages 
to guide us, so we make a mean field approximation for 
the size of the neighbouring gap, namely y = 1. Denot- 
ing the probability of fragmenting a gap into proportions 
a and (I — a) by -P(a), we find the following distribu- 
tional fixed point equation (DFPE) for the probability 
distribution function 4>{x) of gaps x £ [0,oo): 



x = a(l 



(1) 



This convenient notation (exploited below) states that 
the distribution of the variates on the left is equal to that 
on the rights^. The DFPE leads to the Integral Equation 
(IE) for <f>(x), 



4>(x) 



- 1 



P(a) 



da 



(2) 



where we assume that 4>(x) = if x < 0. 

This equation states that the statistical distribution of 
gaps is unchanged by the fragmentation of all the gaps 
incremented in scaled size by one. Note that we neglect 
long-range chronological effects here, of the type appar- 
ent in Figure Q] for the creation of gap 174 which arose 
from the nucleation of island I4 and the fragmentation of 
g a P (53 + 95 + 9i)- We will return to this point below. 

In the aggregation regime, the probability P(a) of frag- 
menting a gap into proportions a and (1— a) is found from 
the steady-state monomer density profile^: 



P(a) 



a Q (l 



B(a + l,a + 1) 
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(3) 



Here B(m,n) = (m — l)!(n — l)!/(m + n — 1)! is the 
Beta function and a reflects the dominant nucleation 
mechanism. For nucleation triggered by deposition of 
monomers, a = i for i = 1,2,3..., whereas for nucle- 
ation resulting from the diffusion of mature monomers 
a = i + l, i = 0,1,2,.. £L. 

In Figure[5]we show the convergence of iterates of equa- 
tion (2), with P(a) given by equation (3). The limit sat- 
isfies the DFPE (1), and so is the form that we wish 
to compare to the scale-inavriant GSD found in the MC 
simulations. 

In Figure |3] we compare the numerically computed 
fixed points of equation @, which we denote by 4> a (x), 
with the GSDs of our MC simulations^ for various crit- 
ical island size i. The comparison is rather good. For 
i = 1,2,3, we see that the observed GSD lies between 
that of the a = i + 1 and a — i distributional fixed point 
solutions. This can be expected since we have found else- 
where that island nucleation is driven by both deposition 
events and purely diffusional fluctuations in monomer 
density^ 7 -. For spontaneous nucleation where i = 0, only 
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FIG. 2. The evolution of gap size distribution 
under iteration of equation ((2]| with i = 1. 
The solid lines are for a = i + 1 in equa- 
tion @, and the broken lines for a — i, where 
the broken lines are shifted along the abscissa 
for clarity. 



the a = i + 1 = 1 model is physically reasonable, since 
there is no possibility of a monomer depositing close to 
a pre-existing critical island of size i in this case. 

It is interesting to ask how the solutions to the DFPE 
equation (J2J) compare to those of the forward-propagated 
fragmentation theory equations, for which the asymp- 
totic behaviours are know n 14 ' 15 . As in [16], it can be 
shown that 4> a {x) obeys the delay-differential equation 
(DDE) 



•2a+l 



d a+1 

dx a+1 



cj> a (x) 



(-l) a+1 a! 
B(a + l,a + l) 



0a(z-l), ( 4 ) 



from which it follows that 



4>a{x) ~ kx C 



as x 



0, 



for some constant k. This is the same small-size 
asymptotic behaviour found in the fragmentation the- 
ory approac h 15 ' 17 . We could not obtain the large-size 
asymptotics for 4> a (x) from the DDE. However, numeri- 
cal analysis of the solutions in Figure [3] show that they 
differ from those obtained by the fragmentation equa- 
tion approach. The reason for this can be traced to the 
derivation of equation ((T|), where not only do we adopt 
a mean field approach for nearest neighbour gap sizes, 
but we also neglect longer-range correlations which are 
expected to be more prominent for larger gaps created 
early in the growth process. Nevertheless, the results in 
Figure [3] show that the solutions capture much of the 
essential physics for the GSDs. 
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We also plot the Generalised Wigner Surmise (GWS) re- 
cently proposed^ as a convenient analytical form for the 
CZDs whose only parameter depends solely on the criti- 
cal island size i. We see that the solution of equation ([5]) 
fits the data at least as well as, and in the case of i = 
much better than, the GWS. Whilst the validity of the 
GWS has been questione d 19 ' 20 , it is a useful benchmark 
for comparisons to MC dat a 21 ' 22 . 




FIG. 3. The GSDs compared to histograms of MC date 
for various critical island size i, taken at nominal coverage 
6 — 100% (9 is deposition rate times elapsed time, and since 
we use point islands, 9 can be much larger than 1). The 
solid curves are the converged solutions to equation (|2j with 
a = i+ 1 with i = 0, 1, 2, 3, and the broken lines are for a — i 
with i — 1,2, 3. 



We turn now to consider the evolution of the capture 
zones in the system. Referring back to Figure [1] we see 
that the capture zone C5 was created by the nucleation 
of island 1$. Prior to this, the zones C3 and C4 were 
larger, so that the creation of C5 can be viewed as the 
fragmentation of part of C3 (the part to the right of is- 
land I 3 ) and part of C4 (to the left of I 4 ). In general 
we do not know how much of the neighbouring capture 
zones to take, nor indeed how large these zones are. How- 
ever, we can again invoke a mean field approximation for 
these nearest neighbour correlations to find the following 
DFPE for a general capture zone c: 




FIG. 4. The CZDs compared to histograms of MC data^ 
for various critical island size i, taken at nominal coverage 
6 — 100%. The solid curves are the solutions to equation ([5]) 
for a — i + 1 with i = 0,1, 2, 3, and the broken lines are for 
a = i with i = 1, 2, 3. 

We can easily quantify the performance of the solutions 
using the moments S m of the distributions. Following 
[23| . from equation © we find the following recursive 
relationship: 



where 



1 \ rn rn 



2 J ^ k\(m-k)\ 



p=0 



p\(m — p)! 



Sp, 
(6) 



A 1 . . - . 

c = 2*- ai + ^2) (1 + Cj. 



(5) 



B„ 



The proportions a\ and 02 are independently drawn from 
P(a) of equation (J3J. An equivalent IE, like that of equa- 
tion ^j, can readily be identified for equation ([5]). 

In Figure |4] we compare the Capture Zone Distribu- 
tions (CZDs) obtained as fixed points of equation ([5]) 
with those from the MC simulations 17 . Again we find 
excellent agreement, particularly for i = and i = 1. 



B(m + a + l,a+l) 
B(a + l,a + l) ' 



and B(m,n) is the Beta function as in Eqn. (|3"T). 
The moments of the GWS are given by 



(< + 1)!' 
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In Table Q] we compare the moments calculated from 
equations ^ and ([7} alongside those taken from our MC 
simulations 17 for i — and i = 1. These confirm the 
superiority of the DFPEs, notably for a = i + 1 implying 
a greater significance for nucleation driven by monomer 
diffusion. 
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1.178 


1.138 


1.134 ±0.001 


1.105 


1.138 


1.098 


1.098 ±0.001 


3 


1.571 


1.439 


1.425 ±0.001 


1.325 


1.439 


1.305 


1.307 ±0.001 


4 


2.313 


1.989 


1.949 ±0.001 


1.708 


1.989 


1.665 


1.666 ±0.001 



TABLE I. Moments of the CZDs from the GWS (equation (0, 
the distributional fixed point equations (equation with 
a — t + 1 or a = z if appropriate), and from the MC 
simulations 17 taken at 8 = 100%. 



In summary, we have presented distributional fixed 
point equations for the nucleation of point islands in one 
dimension. The approach develops a new retrospective 
view of how the inter-island gaps and capture zones have 
developed from the fragmentation of larger entities. That 
the resulting integral equations evolve to fixed points of- 
fers a new perspective on why scale-invariant distribu- 
tions arise in the nucleation process. We have also shown 
that the solutions of our equations compare well to Monte 
Carlo simulation datatf, performing at least as well as the 
recently proposed Generalised Wigner Surmise^, and no- 
tably better for the case of i = 0. 
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